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Abstract. We calculate the rapidity dependence of leading-order medium dilep- 
ton yields resulting from a quark-gluon plasma which has a local time-dependent 
anisotropy in momentum space. We present a phenomenological model which in- 
cludes temporal evolution of the plasma anisotropy parameter, £, and the hard mo- 
mentum scale, fhard- Our model interpolates between 1+1 dimensional collisionally- 
broadened expansion at early times and 1+1 dimensional ideal hydrodynamic ex- 
pansion at late times. Using our model, we find that at LHC energies, forward 
high-energy medium dilepton production would be suppressed by up to a factor 
of 3 if one assumes an isotropization/thermalization time of 2 fm/c. Therefore, 
it may be possible to use forward dilepton yields to experimentally determine 
the time of onset of locally isotropic hydrodynamic expansion of the quark-gluon 
plasma as produced in ultrarelativistic heavy-ion collisions. 



1 Introduction 

Nucleus-nucleus collisions at high energies offer us the opportunity to study experimentally and 
theoretically a possible new state of matter formed by liberated partons known as a quark-gluon 
plasma (QGP). Nevertheless, the main goal is not only to discover this new state of matter but 
also to characterize experimentally many of its properties such as the thermodynamical aspects 
of this new phase. It is very important to know if the matter created after the collision is really 
thermalized. In this direction, an open question is to determine the isotropization and thermal- 
ization time of the matter, Tj SO and r t hcrm, respectively!]] Based on data from the Brookhaven 
National Lab's Relativistic Heavy Ion Collider (RHIC) early studies found that for px < 2 GeV, 
the elliptic flow, V2 , was very well described by models which assumed ideal hydrodynamic be- 
havior starting at very early times r ~ 0.6 fm/c [l|2|3j . These fits of v-i indicated that the 
matter could be modeled as a nearly-perfect fluid and hence implied fast thermalization of the 
matter created at RHIC energies. However, recent results from conformal relativistic viscous 
hydrodynamics [4 have shown that these initial estimates for the isotropization/thermalization 
time of the plasma are not completely reliable due to poor knowledge of the proper initial 
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Hereafter, for simplicity, will assume that these two time scales are the same so that Ttherm — Tiso- 
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conditions (CGC versus Glauber), details of plasma hadronization such as the choice of the 
proper freezeout time and the subsequent hadronic cascade, etc. Now it seems that isotropiza- 
tion times up to Tj SO ~ 2 fm/c are not ruled out by RHIC data. Therefore, additional theoretical 
and experimental input are necessary to further constrain this time. 

One way to increase our understanding of the pre-equilibrium stage of QGP evolution is 
by studying independent observables which are sensitive to early-time dynamics. One good 
candidate is high-energy dileptonsB Due to their large mean free path, lepton pairs can leave 
the strongly interacting region carrying information about the earliest times after the nuclear 
collision. In this work, we calculate the rapidity dependence of dilepton pair production from a 
QGP which has a time-dependent anisotropy in momentum space. 

Most previous phenomcnological treatments of dilepton production have assumed that the 
plasma is thermalized rapidly with t iso on the same order as the parton formation time, tq. How- 
ever, due to the rapid 1+1 longitudinal expansion of the plasma at early-times, this assumption 
seems to be rather strong because it ignores the momentum-space anisotropy developed along 
the beam axis. Anisotropies in momentum-space are intimately related with Weibel instabilities. 
Recently, it has been shown that the presence of chromomagnetic Weibel instabilities developed 
at early times of the collision may affect the value of Ti so [5I6I7I8I9] , However, it is still not clear 
how much Ti SO is affected by the presence of these instabilities at phenomenologically accessible 
collision energies. 

Due to the lack of a precise physical picture of the pre-equilibrium QGP, one can instead im- 
plement simple models which capture the essence of the complex dynamics at early-times of the 
collision. This will be our approach in the present work. We propose a simple phenomcnological 
model for the time-dependence of the hard momentum scale, Phard, and the plasma anisotropy 
parameter, £ = (j>j,) /2{p 2 L ) — 1. The new aspect in this work compared to our recent proposals 
for the pre-equilibrium stage of the QGP |10|11|12] , is the inclusion of the rapidity dependence 
of the quark and anti-quark distribution functions. We use fits to experimental data which are 
available from AGS through RHIC data to constrain the rapidity dependence of the parton 
distribution functions. Since this rapidity dependence is not perfectly flat we implicitly include 
effects of the breaking of longitudinal boost invariance of the system. 

To make our final phenomcnological predictions for dilepton yields, we introduce three 
different parameters: (1) the parton formation time, To, which is time at which the partons 
created during the initial hard collision of the nuclei become decoherent and approximately on- 
shell0 (2) Tj SO which is the proper time when the system starts to undergo ideal hydrodynamical 
expansion; and (3) 7 which sets the sharpness of the transition to hydrodynamical behaviour. 
We assume that for times of the order of To but smaller than Ti so , the system expands like a 
collisionally-broadened plasma and smoothly changes to hydrodynamical expansion for times 
long compared with t iso . In this work we extend previous calculations [10111112"] considering 
the impact of the momentum-space anisotropies on the full rapidity dependence of medium 
dilepton production. We find that at LHC energies, dilepton yields are suppressed by a factor 
of 3 around y ~ 9 if one chooses Tj SO =2 fm/c. 

This manuscript is organized as follows: In Sec. [2] we calculate the dilepton production 
rate at leading order using an anisotropic phase space distribution. In Sec. [3] we present the 
interpolating model from collisionally-broadened to ideal hydrodynamical expansion including 
the rapidity dependence of the quark and anti-quark distribution functions. In Sec. [4] we present 
dilepton yields as a function of the rapidity for different values of Ti SQ . Finally, we present our 
conclusions and give an outlook in the Sec.[5J 



2 We mean by high-energy dileptons, lepton pairs with pair transverse momentum (pr) or invariant 
mass (M) greater than 1 GeV. 

3 This time can be estimated from the nuclear saturation scale, i.e., to ~ QJ 1 . For LHC energies, 
Q s ~ 2 GeV implying tq ~ 0.1 fm/c 
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2 Dilepton rate 

From relativistic kinetic theory, the dilepton production rate dN l+l jd 4 Xd 4 P = dR l+l /d 4 P 
(i.e. the number of dileptons produced per space-time volume and four dimensional momentum- 
space volume) at leading order in the electromagnetic coupling, a, is given by |14|15|16] : 



dR l+l _ f d 3 P! d 3 p 2 



d*P~-J (2^)3 ^fMf*fa) v rt a W~ 6U) ( p -P* (!) 

where f q< q is the phase space distribution function of the medium quarks (anti-quarks), v qq - is 
the relative velocity between quark and anti-quark and a qq is the total cross section. Since 
we will be considering high-energy dilepton pairs with center-of-mass energies much greater 
than the dilepton mass we can ignore the finite dilepton mass corrections and use simply 
a qq = 4ira 2 /3M 2 . In addition, we assume that the distribution function of quarks and anti- 
quarks is the same, f q = f q . 

Following Ref. [5], we will consider the case that the anisotropic quark/anti-quark phase- 
space distributions are azimuthally symmetric in momentum space about a direction specified 
by n and can be obtained from an arbitrary isotropic phase space distribution by squeezing 
(£ > 0) or stretching (£ < 0) this isotropic distribution function along n, i.e., 



/,,*(p, £,Phard) = /^U/p 2 +£(p-fi) 2 ,Phard) , (2) 

where phard is the hard momentum scale, n is the direction of the anisotropy and £ > is 
a parameter that reflects the strength and type of anisotropy. In the isotropic case, Phard is 



identified with the temperature of the system and £ = 0. Using the anisotropic distribution 
given by Eq. © into Eq. (TT]) we obtainQ 

dR l+r 5a 2 f 1 „ „ , f a - d V i 



^d(cos0 pi ) J -j=pi f q ^yp?(l +C(r)cos 2 e» Pl ),p hard (T,7 / ) 



with 



fq- x /(E-p 1 )2 + £(r)(p 1 cos0p 1 -Pcos0 P )2,p hard (r,^) , (3) 



X = AP 2 p\ sm 2 9p sin 2 ^ - {2 Pl {E - PcosO P cos6» Pl ) - M 2 ) 2 
M 2 



a± = 



2{E - Pcos(9 P ±e pi )) 



For phenomenological predictions of expected dilepton yields, we model the space-time de- 
pendence of Phard and £. In a realistic non-boost-invariant system phard will depend of the 
rapidity on the quark and anti-quark and hence the rate itself will depend not only on the 
difference of y — rj but also on rj itself. In the next sections we will give details of how this de- 
pendence is introduced. Once obtained we can calculate the final rapidity spectra by integrating 
over phase space: 



^L— = *R 2 T j dM 2 d 2 vr J J ^qrrdrdr, , (4) 



where Rj- = 1.2 A 1 / 3 fm is the radius of the nucleus in the transverse plane. This expression 
is evaluated in the center-of-mass (CM) frame while the differential dilepton rate is calculated 
in the local rest frame (LR) of the emitting region. Then, the dilepton pair energy has to be 

understood as Elr = Pt cosh (y — 77) in the differential dilepton rate dR l 1 /d 4 P. Substituting 



4 Details of the calculation are presented in Ref. |12] , 
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Eq. ((3]) into Eq. (j4|) , we obtain the dilepton spectrum as a function of the rapidity including 
the effect of a time-dependent momentum anisotropy. 

One can be worried if either transverse expansion or mixed/hadronic phases in Eq. Q will 
affect the production of high-energy dileptons presented here. Fortunately, in the kinematic 
region studied here this is not the case and these effects turn out to be negligible (1-2% effect) 
compared with the longitudinal expansion [10] . 



3 Space Time-Model 

In this section, we extend previous models proposed for the pre-equilibrium stage of the QGP 
|10|11|12] by including a phcnomenologically realistic rapidity dependence of the parton dis- 
tributions. Before presenting the details of the model, we remind the reader of some general 
relations. 

Firstly, the plasma anisotropy parameter is related with the average longitudinal and trans- 
verse momentum of the hard particles through the relation: 

-«> (5) 



s 2<pD ■ 

We can obtain two limiting cases from the last expression. For an isotropic plasma, we have 
that £=0, i.e., {p\) = 2(p|_). Another possibility is that the partons initially undergo Id free 
streaming expansion, where the partons expand freely along the longitudinal axis. Using the 
free streaming distribution function, it is possible to show that the transverse and longitudinal 
momentum scales as |9I12I13| : 

(p 2 T } f . s . <x2T 2 , (6a) 

(Pi)f.,«^. (6b) 

Inserting these expressions into Eq. ([5]), one obtains £/. s .(t) = t 2 /tq — 1. The free streaming 
result for £ is the upper-bound on possible momentum-space anisotropies developed during Id 
expansion by causality. Modifications to the Eq. ^ result from the different kinds of interactions 
between the partons, as it is discussed below. In this work, for simplicity, we will not study 
explicitly the possibility of Id free streaming expansion since, in reality, this is a rather extreme 
assumption which requires that the partons do not interact at all. Below we will give details of a 
more realistic model of the time evolution of the average transverse and longitudinal pressures. 
However, note that the relation given in Eq. is completely general and can be applied in all 
cases. 

Secondly, for a given anisotropic phase space distribution of the form specified in Eq. ([2| , 
the local energy density can be factorized as: 

£0hard, = J 7^X3 P fiso ( VP 2 + I (P ' n 2 ) , Phard ) , (7) 
= So(Phard)ft(£) ) 

where £q is the initial local energy density deposited in the medium at To and 

R(0 4(_i_ + «). (8) 

We mention that phard depends implicitly on proper time and space-time rapidity and that 
lim^o TZ{() = 1 and lim^^ = 1/ vt 
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3.1 Momentum-space broadening and plasma instabilities effect 

The ratio between the average longitudinal and transverse momentum needed to compute £ 
using Eq. ([5]) is modified from the free streaming case if collisions between the partons are 
taken into account. To include all relevant N <-> N Boltzmann collision terms plus mean field 
interactions (Vlasov term) is a very difficult task. As a first approximation, one can consider 
the effect of elastic collisions in the broadening of the longitudinal momentum of the particles. 
This was the approach in the original version of the bottom-up scenario [13] . In the first stage 

3/2 

of this scenario, 1 <C Q s t <C a s , initial hard gluons have typical momentum of order Q s 
and occupation number of order l/a s . Due to the fact that the system is expanding at the 
speed of light in the longitudinal direction, the density of hard gluons decreases with time, 
N g ~ Qs/(a s Q s T). If there were no interactions this expansion would be equivalent to 1+1 free 
streaming and the longitudinal momentum pi, would scale like 1/r. However, once elastic 2 <-» 2 
collisions of hard gluons are taken into account, the ratio between the longitudinal momentum 
Pl and the typical transverse momentum of a hard particle px decreases as |13j : 

W .r)^. (9) 

Assuming, as before, isotropy at the formation time, tq — Qj 1 , this implies that for a collisionally- 
broadened expansion, £(t) = (t/t ) 2 / 3 — 1. 

Note that, as obtained in Ref. 13J, when collisions are included, it is implicitly assumed that 
the elastic cross-section is screened at long distances by an isotropic real- valued Debye mass. 
This is not guaranteed in an anisotropic plasma since the Debye mass can be become complex 
due to the chromo-Weibel instability. In general, any anisotropic distribution causes the presence 
of negative eigenvalues in the self-energy of soft gluons [5I6I7I8I9] , Such negative eigenvalues 
indicate instabilities, which result in exponential growth of chromo-electric and magnetic fields, 
E a and B a , respectively. These fields deflect the particles and how much deflection occurs will 
depend on the amplitude and domain size of the induced chromofields. Currently, the precise 
parametric relation between the plasma anisotropy parameter and the amplitude and domain 
size of the chromofields is not known from first principles. 

If one would like to include the momentum-space broadening and the impact of the plasma 
instabilities, phenomenologically this can be achieved by generalizing the temporal dependence 
of £(t) to: 

S 



The exponent S contains the physical information about the particular type of momentum-space 
broadening which occurs due to plasma interactions. Two limiting cases for this exponent are 
the ideal hydrodynamic and free streaming expansion. In the 1+1 hydrodynamical limit, S = 
and then £ — > 0. For 5=2, one recovers the 1+1 dimensional free streaming case, £ — > £f. s . = 
(t/to) 2 — 1. For general S between these limits one obtains the proper-time dependence of 
the energy density and temperature by substituting (|10[) into the general expression for the 
factorized energy density j7|) to obtain £(t,5) = £o(j?hard) 7?-(£(t, <5)). The inclusion of the 
rapidity dependence of the parton distribution functions is discussed in Sec. 13.21 

Different values of 6 arise dynamically from the different processes contributing to parton 
isotropization. For example, elastic collisional-broadening results in Eq. ((9J and hence S = 2/3. 
Recently, some authors have considered the values of 6 resulting from processes associated with 
the chromo-Weibel instability presented at the earliest times after the initial nuclear impact 
[17118119] : 



where 



Ref. Q2] , 

1 Ref. [E] , (12) 

2 Nielsen-Olesen limit, Ref. [19] . 
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Fig. 1. Sketch of the time dependence the plasma anisotropy indicating the various time scales and 
processes taking place. Here tmfp is the mean time between elastic collisions (mean-free time) and 
^instability is the time at which plasma-instability induced soft modes have grown large enough to affect 
hard particle dynamics. 



Summarizing, the coefficient <5 in Eq. (|10p takes on the following values 

2 Free streaming expansion , 
2/3 Collisional-Broadening, Ref.[T3] , 
l/2Ref.[Tj], 

1/4 Ref.pl , 1 J 

1/6 Nielscn-Olcsen limit, Ref.Ql] , 
Hydrodynamic expansion . 

In Fig. [1] we sketch the time-dependence of the plasma anisotropy parameter indicating the 
time scales at which the various processes become important. At times shorter than the mean 
time between successive elastic scatterings, tmfp, the system will undergo 1+1 dimensional 
free streaming with 6 = 2. For times long compared to tmfp but short compared to Testability 
the plasma anisotropy will grow with the collisionally- broadened exponent of 6 = 2/3. Here 
^instability is the time at which instability-induced soft gauge fields begin to influence the hard- 
particles' motion. When Ti nsta biiity < t < Tj SO the plasma anisotropy grows with the slower 
exponent of S = 1/6 ... 1/2 due to the bending of particle trajectories in the induced soft-field 
background. At times large compared to Testability inelastic processes are expected to drive 
the system back to isotropy [13]. We note here that for small £ and realistic couplings it has 
been shown [20] that one cannot ignore the effect of collisional-broadening of the distribution 
functions and that this may completely eliminate unstable modes from the spectrum. 

Based on such a sketch, one could try to construct a detailed model which includes all of 
the various time scales and study the dependence of the process under consideration on each. 
However, there are theoretical uncertainties in each of these time scales and their dependences 
on experimental conditions. We choose to use a simpler approach in which we will construct a 
phcnomcnological model which smoothly interpolates the coefficient 5 from the Id collisionally- 
broadened expansion to Id hydrodynamical expansion, i.e., 2/3 > 6 > 0. 

In the model we introduce a transition width, 7 , which governs the smoothness of the 
transition from the initial value of<5 = 2/3to(5 = OatT~ Ti SO - The collisionally-broadened 
interpolating model provides us a realistic estimate of the effect of plasma anisotropies. Note 
that by using such a smooth interpolation one can achieve a reasonable phenomenological de- 
scription of the transition from non-equilibrium to equilibrium dynamics which should hopefully 
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capture the essence of the physics. In the next section we will give mathematical definitions for 
the model. 



3.2 Interpolating model for collisionally broadened expansion 

In order to construct an interpolating model between collisionally-broadened and hydrodynam- 
ical expansion, we introduce the smeared step function: 



A(t, T iso ,7) = i ^tanh 



7(7 



1 



(14) 



where 7 -1 sets the width of the transition between non-equilibrium and hydrodynamical evo- 
lution in units of Tj S0 . In the limit when r <C t iso , we have A — * and when t Tj SO we have 
A — * 1. Physically, the energy density £ should be continuous as we change from the initial 
non-equilibrium value of <5 to the final isotropic 5 — value appropriate for ideal hydrodynamic 
expansion. Once the energy density is specified, this gives us the time dependence of the hard 
momentum scale. We find that for general 8 this can be accomplished with the following model: 



with 7£(£) defined in Eq. 



e(r,<5) = (r/ro) <5(1 - A(T)) -l J 

Phard(T, J? )=ToW 1 / 3 (r) F(r?), 
I and for fixed final multiplicity we have: 



(15a) 
(15b) 
(15c) 



- 1 



U(t) = [^((T iso /r )^ 
U(t)=U(t)/U(t+), 



U{t^ )= lim W(r) = Tl 



3A(t)/4 



7~iso 

r 



l-A-(l-A(r))/2 



((Tiso/ro) 15 



3/4 



7"0 



(16a) 
(16b) 



and (5 = 2/3 for the case of Id collisionally broadened expansion interpolating to Id ideal 
hydrodynamic expansion. 

Recently, models for the pre-equilibrium stage in the presence of momentum-space anisotropics 
have been applied in the central rapidity region of high-energy dileptons [10111112] . Here we 
extend these previous analyses through the inclusion of rapidity dependence of the parton hard 
momentum scale by using a phenomenologically constrained "rapidity profile function," F(rj). 
Our main goal with this modification is to explore the phenomenological consequences in the 
forward rapidity region where the effect of early-time anisotropies are expected to be maxi- 
mal. We are not attempting to describe the physics of the forward rapidity region from first 
principles^ instead, we implement a Gaussian fit profile for the rapidity dependence which suc- 
cessfully describes experimentally observed pion rapidity spectra from AGS to RHIC energies 
|25l26l27l28l29j and use this to extrapolate to LHC energies: 



(17) 



with 



3(l-cf) 

where c s is the sound velocity and m p is the proton mass. 



(18) 



5 Some proposals have been mentioned in the literature, see Ref. 121122123124] 
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Fig. 2. Temporal evolution of the plasma anisotropy parameter using our collisionally-broadened 
interpolating model for four different isotropization times n so £ {1, 4, 6, 18} to. The transition width is 
taken to be 7 = 2. To convert to physical scales use to ~ 0.1 fm/c for LHC energies. 



Note that once the rapidity dependence in the parton momentum distribution functions is 
implemented boost invariance along the longitudinal axis breaks down. This procedure leads 
to a violation of the conservation laws expressed by hydrodynamics unless a finite baryon 
chemical potential is introduced [30] . It is possible to demonstrate that for a longitudinal scaling 
expansion, dV(T, fj.) /di]=0, where V is the pressure and T is the temperature. This condition 
is equivalent to(£| 

dT da 

s ^ +n i=^ (19) 

where s and n denote the entropy density and particle number density, respectively. In the 
present context don't have to worry about the presence of finite chemical potentials since 
we are considering high-energy dilepton production and E/T ^> /x/T is satisfied. Therefore, 
in the differential dilepton rate dR/d^P, Eq.JT]), the product of the distribution functions 

/g(Pg,+M) /«(P?,-M) ~ fq(Pq) fq(Pq)- 

In Fig.O the temporal evolution of the anisotropy parameter £(r) is plotted using Eq. (|15a[) . 
In Fig.[3j we show the time and rapidity dependence of Phard(T, v) (right and left panel, respec- 
tively) using Eq. (|15c|) . 



4 Results 

In this section, we will present our predicted dilepton yields as a function of the rapidity from 
a Pb-Pb collision at LHC full beam energy, y / SNN= 5.5 TeV. At this center-of-mass energy 
we use t = 0.088 fm/c, T = 845 MeV, R T = 7.1 fm and the critical temperature T c =160 MeV. 
The kinematic cuts in the transverse momentum and invariant mass of the dilepton yields are 
indicated in the corresponding results. Also, we use (? s = 1/3 and m p — 0.938 GeV in Eq. (jT5J). 

Before presenting our results we first explain the numerical procedure used for our calcula- 
tions. Because the differential dilepton rate dR l 1 /d 4 P given in Eq. ([3]) is independent of the 
assumed space-time model, we first evaluate it numerically using double-exponential integra- 
tion with a target precision of 10 -9 . The result for the rate was then tabulated on a uniformly- 
spaced 4-dimensional grid in M, px, y, and £: M/phard,_PT/Phard € {0.1, 20}, y G { — 10, 10} and 



Partial derivatives with respect to r\ are performed at constant r. 
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Fig. 3. Temporal evolution (left panel) and rapidity dependence (right panel) of the hard momentum 
scale, Phard using our fixed final multiplicity collisionally-broadened interpolating model for the hard 
momentum scale for four different isotropization times Tiso G {1, 4, 8, 16} to. The transition width is 
taken to be 7 = 2. To convert to physical units use to ~ 0.1 fm/c for LHC energies. For the rapidity 
dependence of Phard (right panel) we used a constant value of r ~ 0.1 fm/c and the width ~ 8. 



£ € {0,5}. This table was then used to build a four-dimensional interpolating function which 
was valid at continuous values of these four variables. We then boost this rate from the local 
reference frame to center-of-mass frame and evaluate the remaining integrations over space-time 
(r and 77), transverse momentum and invariant mass appearing in Eq. (j4]) using quasi-Montc 
Carlo integration with r € {to,t/}, 77 6 {—10,10} and, depending on the case, restrict the 
integration to any cuts specified in M or px- 

Our final integration time, r/, is set by solving numerically for the point in time at which the 
temperature in our interpolating model is equal to the critical temperature, i.e. Pha.rd(i~f , rj) = 
Tc- We will assume that when the system reaches Tc, all medium emission stops. Note that 
due to the fact that phard depends on the parton rapidity, the plasma lifetime now depends on 
which rapidity slice you are in, with higher rapidities having a shorter lifetime due to their lower 
initial "temperature". We are not taking into account the emission from the mixed/hadronic 
phase at late times since the kinematic regime we study (high M and pr) is dominated by 
early-time high-energy dilepton emission [10|16j . 

We show our predicted dilepton spectrum as a function of the pair rapidity, y, for LHC 
energies using our model described by Eqs. (fl~5f in Fig.[4j From this, we see that for LHC energies 
there is a suppression when we vary the isotropization time from To to 2 fm/c. This suppression 
can be explained qualitatively by two mechanisms. The first one, the anisotropic nature of 
the distribution function as a consequence of the rapid expansion implies that dileptons with 
larger values of longitudinal momentum are reduced compared with the case of an isotropic 
distribution function. The suppression will depend on the maximum amount of momentum- 
space anisotropy achieved at early times and also on the time dependence of the anisotropy 
parameter £; in this work, we consider a realistic scenario for a collisionally-broadened plasma. 
The other source of rapidity dependence of the final dilepton spectra is related to the fact that 
the hard momentum scale ( "temperature" ) depends explicitly on the rapidity 77, even in the 
case of instantaneous thermalization. To generate Fig. 2] we have applied a cut M > 2 GeV 
and Pt > 100 MeV. As can be seen from this figure a isotropization time of Tj SO = 2 fm/c 
results in fewer dileptons as compared to "instantaneous" isotropization Tj so = 0.088 fm/c. 
This suppression is enhanced at forward rapidities. 

In order to quantify the effect of the pre-equilibrium emission we define the "dilepton mod- 
ification" factor as the ratio of the dilepton yield obtained with an isotropization time of Ti so 
to that obtained from an instantaneously thermalized plasma undergoing only 1+1 hydrody- 
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Fig. 4. Fixed final multiplicity condition collisionally-broadened interpolating model dilepton yields 
as a function of rapidity in Pb-Pb collisions at LHC, with a cut M > 2 GeV and Pt > 100 MeV. For 
medium dileptons we use 7 = 2 and Ti so to be either 0.088 or 2 fm/c for LHC energies. 



namical expansion, ie. t- iso = tq: 

. (^^L\ I (^z^lA\ . (20) 

This ratio measures how large the effect of early-time momentum anisotropies are on medium 
dilepton production. In the case of instantaneous isotropization, <P(Ti so ) is unity, and for t; so > To 
any deviation from unity indicates a modification of medium dilepton production due to pre- 
equilibrium emissions. 

In Fig. [5] we show our prediction for the rapidity dependence of the high-energy dilepton 
modification factor, <£(ri so ), for three different assumed plasma isotropization times, t; so G { 
0.1,1,2} fm/c. To generate this figure we have applied a cut M > 2 GeV and Pt > 100 MeV. 
As can be seen from this figure a isotropization time of t; so = 2 fm/c results in fewer dileptons 
as compared to "instantaneous" isotropization t iso = 0.088 fm/c. This suppression is enhanced 
at forward rapidities and reaches a maximum suppression of a factor of 3 at extremely forward 
rapidities. 

Using the dilepton modification factor as our criterion we find that for our collisionally- 
broadened interpolating model with fixed final multiplicity, the dilepton yields as a function of 
the rapidity at Ti SO = 2 fm/c can be suppressed up to ~ 20% for 0< y < 4. The suppression of 
dilepton yields is more dramatic at rapidity values around y ~ 9 and can be as large as a factor 
of 3. With sufficiently accurate experimental results this could give an experimental method for 
determining the isotropization time of a quark gluon plasma as formed in an ultrarelativistic 
nuclear collision. 



5 Conclusions 

In this work we have introduced a phcnomcnological model that takes into account early-time 
momentum-space anisotropies in the rapidity dependence of high-energy dilepton production. 
To do this we have modeled the temporal evolution of the plasma anisotropy parameter £ and 
the hard momentum scale Phard- To study the dilepton production rapidity dependence, we 
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Fig. 5. Predicted dilepton modification factor, <P(ri so ), for three different assumed plasma isotropization 
times, Ti so e { 0.1,1,2} fm/c. Cuts are the same as in Fig. [4] 

have parametrized the rapidity dependence of phard using a Gaussian profile which is consistent 
with experimental observations of final pion spectra from AGS through RHIC energies. 

We have applied the proposed model to study high-energy dilepton yields as a function 
of the pair rapidity and find that this observable is sensitive to the chosen value of t- iso . This 
suppression can be explained as a consequence of the combined effect of the anisotropy in 
momentum-space achieved at early-times due to expansion and the rapidity dependence of the 
hard momentum scale which explicitly breaks longitudinal boost invariance. We find that with 
the resulting dilepton modification factor, <P(Ti so =2 fm/c), shows suppressed dilepton yields in 
the forward rapidity region which can be up to 20% for < y < 4 and up to a factor of 3 at 
y ~ 9. The amplitude of the suppression of ^(r iso ) could help us to experimentally constrain 
Ti so given sufficiently precise data in the forthcoming LHC experiments. In this way forward 
dileptons would provide a way to determine the plasma isotropization time experimentally. 

An uncertainty of our treatment comes from our implicit assumption of chemical equilib- 
rium. If the system is not in chemical equilibrium (too many gluons and/or too few quarks) 
early time quark chemical potentials, or fugacities, will affect the production of lepton pairs 
[15116) . However, to leading order, the quark and gluon fugacities will cancel between numerator 
and denominator in the dilepton suppression factor, $(Ti so ) [E]. We, therefore, expect that to 
good approximation one can factorize the effects of momentum space anisotropics and chemical 
non-equilibrium. 

We note in closing that the interpolating model presented here can be applied to other ob- 
servables than dilepton yields. Indeed, with this model it is possible to assess the phenomeno- 
logical consequences of momentum-space anisotropies on other observables which are sensitive 
to early-time stages of the QGP, e.g. photon production, heavy-quark transport, jet-medium 
induced electromagnetic and gluonic radiation, etc. 
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